Multifaceted immune dysregulation characterizes individuals at-risk for rheumatoid arthritis.

Molecular markers of autoimmunity, such as antibodies to citrullinated protein antigens (ACPA), are detectable prior to inflammatory arthritis (IA) in rheumatoid arthritis (RA) and may define a state that is 'at-risk' for future RA. Here we present a cross-sectional comparative analysis among three groups that include ACPA positive individuals without IA (At-Risk), ACPA negative individuals and individuals with early, ACPA positive clinical RA (Early RA). Differential methylation analysis among the groups identifies non-specific dysregulation in peripheral B, memory and naïve T cells in At-Risk participants, with more specific immunological pathway abnormalities in Early RA. Tetramer studies show increased abundance of T cells recognizing citrullinated (cit) epitopes in At-Risk participants, including expansion of T cells reactive to citrullinated cartilage intermediate layer protein I (cit-CILP); these T cells have Th1, Th17, and T stem cell memory-like phenotypes. Antibody-antigen array analyses show that antibodies targeting cit-clusterin, cit-fibrinogen and cit-histone H4 are elevated in At-Risk and Early RA participants, with the highest levels of antibodies detected in those with Early RA. These findings indicate that an ACPA positive at-risk state is associated with multifaceted immune dysregulation that may represent a potential opportunity for targeted intervention.

Investigating iron intake in risk of progression from islet autoimmunity to type 1 diabetes: The diabetes autoimmunity study in the young.

Background: Studies of the role of iron in the risk of type 1 diabetes (T1D) have been inconsistent. Given that iron generates reactive oxygen radicals, which can lead to oxidative damage and apoptosis in the beta cells of the pancreas, we examined whether iron intake was associated with the risk of progressing to T1D in individuals with islet autoimmunity (IA), the pre-clinical phase of T1D. Methods: DAISY is a prospective cohort following 2,547 children at increased risk for IA and progression to T1D. IA is defined as at least two consecutive serum samples positive for at least one autoantibody (insulin, GAD, IA-2, or ZnT8). We measured dietary intake at the time of IA seroconversion in 175 children with IA, and of these, 64 progressed to T1D. We used Cox regression to examine the association between energy-adjusted iron intake and progression to T1D, adjusting for HLA-DR3/4 genotype, race/ethnicity, age at seroconversion, presence of multiple autoantibodies at seroconversion, and multiple vitamin use. In addition, we tested whether this association was modified by vitamin C or calcium intake. Results: In children with IA, high iron intake (as defined as above the 75th percentile, > 20.3 mg/day) was associated with decreased risk of progression to T1D compared to moderate iron intake (as defined by the middle 25-75th percentiles, 12.7-20.3 mg/day) (adjusted hazard ratio (HR): 0.35; 95% confidence interval (CI): 0.15, 0.79). The association between iron intake and T1D was not modified by vitamin C nor calcium intake. In a sensitivity analysis, the removal of six children who had been diagnosed with celiac disease prior to IA seroconversion did not affect this association. Conclusion: Higher iron intake at the time of IA seroconversion is associated with a lower risk of progression to T1D, independent of multivitamin supplement use. Further research that includes plasma biomarkers of iron status is needed to investigate the relationship between iron and the risk of T1D.

Genome-wide analysis of oxylipins and oxylipin profiles in a pediatric population.

Background: Oxylipins are inflammatory biomarkers derived from omega-3 and-6 fatty acids implicated in inflammatory diseases but have not been studied in a genome-wide association study (GWAS). The aim of this study was to identify genetic loci associated with oxylipins and oxylipin profiles to identify biologic pathways and therapeutic targets for oxylipins. Methods: We conducted a GWAS of plasma oxylipins in 316 participants in the Diabetes Autoimmunity Study in the Young (DAISY). DNA samples were genotyped using the TEDDY-T1D Exome array, and additional variants were imputed using the Trans-Omics for Precision Medicine (TOPMed) multi-ancestry reference panel. Principal components analysis of 36 plasma oxylipins was used to capture oxylipin profiles. PC1 represented linoleic acid (LA)- and alpha-linolenic acid (ALA)-related oxylipins, and PC2 represented arachidonic acid (ARA)-related oxylipins. Oxylipin PC1, PC2, and the top five loading oxylipins from each PC were used as outcomes in the GWAS (genome-wide significance: p < 5×10-8). Results: The SNP rs143070873 was associated with (p < 5×10-8) the LA-related oxylipin 9-HODE, and rs6444933 (downstream of CLDN11) was associated with the LA-related oxylipin 13 S-HODE. A locus between MIR1302-7 and LOC100131146, rs10118380 and an intronic variant in TRPM3 were associated with the ARA-related oxylipin 11-HETE. These loci are involved in inflammatory signaling cascades and interact with PLA2, an initial step to oxylipin biosynthesis. Conclusion: Genetic loci involved in inflammation and oxylipin metabolism are associated with oxylipin levels.

Propensity scores as a novel method to guide sample allocation and minimize batch effects during the design of high throughput experiments.

BACKGROUND: We developed a novel approach to minimize batch effects when assigning samples to batches. Our algorithm selects a batch allocation, among all possible ways of assigning samples to batches, that minimizes differences in average propensity score between batches. This strategy was compared to randomization and stratified randomization in a case-control study (30 per group) with a covariate (case vs control, represented as ß1, set to be null) and two biologically relevant confounding variables (age, represented as ß2, and hemoglobin A1c (HbA1c), represented as ß3). Gene expression values were obtained from a publicly available dataset of expression data obtained from pancreas islet cells. Batch effects were simulated as twice the median biological variation across the gene expression dataset and were added to the publicly available dataset to simulate a batch effect condition. Bias was calculated as the absolute difference between observed betas under the batch allocation strategies and the true beta (no batch effects). Bias was also evaluated after adjustment for batch effects using ComBat as well as a linear regression model. In order to understand performance of our optimal allocation strategy under the alternative hypothesis, we also evaluated bias at a single gene associated with both age and HbA1c levels in the 'true' dataset (CAPN13 gene). RESULTS: Pre-batch correction, under the null hypothesis (ß1), maximum absolute bias and root mean square (RMS) of maximum absolute bias, were minimized using the optimal allocation strategy. Under the alternative hypothesis (ß2 and ß3 for the CAPN13 gene), maximum absolute bias and RMS of maximum absolute bias were also consistently lower using the optimal allocation strategy. ComBat and the regression batch adjustment methods performed well as the bias estimates moved towards the true values in all conditions under both the null and alternative hypotheses. Although the differences between methods were less pronounced following batch correction, estimates of bias (average and RMS) were consistently lower using the optimal allocation strategy under both the null and alternative hypotheses. CONCLUSIONS: Our algorithm provides an extremely flexible and effective method for assigning samples to batches by exploiting knowledge of covariates prior to sample allocation.

COVID-19 Mortality in the Colorado Center for Personalized Medicine Biobank.

Over 6.37 million people have died from COVID-19 worldwide, but factors influencing COVID-19-related mortality remain understudied. We aimed to describe and identify risk factors for COVID-19 mortality in the Colorado Center for Personalized Medicine (CCPM) Biobank using integrated data sources, including Electronic Health Records (EHRs). We calculated cause-specific mortality and case-fatality rates for COVID-19 and common pre-existing health conditions defined by diagnostic phecodes and encounters in EHRs. We performed multivariable logistic regression analyses of the association between each pre-existing condition and COVID-19 mortality. Of the 155,859 Biobank participants enrolled as of July 2022, 20,797 had been diagnosed with COVID-19. Of 5334 Biobank participants who had died, 190 were attributed to COVID-19. The case-fatality rate was 0.91% and the COVID-19 mortality rate was 122 per 100,000 persons. The odds of dying from COVID-19 were significantly increased among older men, and those with 14 of the 61 pre-existing conditions tested, including hypertensive chronic kidney disease (OR: 10.14, 95% CI: 5.48, 19.16) and type 2 diabetes with renal manifestations (OR: 5.59, 95% CI: 3.42, 8.97). Male patients who are older and have pre-existing kidney diseases may be at higher risk for death from COVID-19 and may require special care.

An Oxylipin-Related Nutrient Pattern and Risk of Type 1 Diabetes in the Diabetes Autoimmunity Study in the Young (DAISY).

Oxylipins, pro-inflammatory and pro-resolving lipid mediators, are associated with the risk of type 1 diabetes (T1D) and may be influenced by diet. This study aimed to develop a nutrient pattern related to oxylipin profiles and test their associations with the risk of T1D among youth. The nutrient patterns were developed with a reduced rank regression in a nested case-control study (n = 335) within the Diabetes Autoimmunity Study in the Young (DAISY), a longitudinal cohort of children at risk of T1D. The oxylipin profiles (adjusted for genetic predictors) were the response variables. The nutrient patterns were tested in the case-control study (n = 69 T1D cases, 69 controls), then validated in the DAISY cohort using a joint Cox proportional hazards model (n = 1933, including 81 T1D cases). The first nutrient pattern (NP1) was characterized by low beta cryptoxanthin, flavanone, vitamin C, total sugars and iron, and high lycopene, anthocyanidins, linoleic acid and sodium. After adjusting for T1D family history, the HLA genotype, sex and race/ethnicity, NP1 was associated with a lower risk of T1D in the nested case-control study (OR: 0.44, p = 0.0126). NP1 was not associated with the risk of T1D (HR: 0.54, p-value = 0.1829) in the full DAISY cohort. Future studies are needed to confirm the nested case-control findings and investigate the modifiable factors for oxylipins.

Hippocampal RNA sequencing in mice selectively bred for high and low activity.

High and Low Activity strains of mice were bidirectionally selected for differences in open-field activity (DeFries et al., 1978, Behavior Genetics, 8: 3-13) and subsequently inbred to use as a genetic model for studying anxiety-like behaviors (Booher et al., 2021, Genes, Brain and Behavior, 20: e12730). Hippocampal RNA-sequencing of the High and Low Activity mice identified 3901 differentially expressed protein-coding genes, with both sex-dependent and sex-independent effects. Functional enrichment analysis (PANTHER) highlighted 15 gene ontology terms, which allowed us to create a narrow list of 264 top candidate genes. Of the top candidate genes, 46 encoded four Complexes (I, II, IV and V) and two electron carriers (cytochrome c and ubiquinone) of the mitochondrial oxidative phosphorylation process. The most striking results were in the female high anxiety, Low Activity mice, where 39/46 genes relating to oxidative phosphorylation were upregulated. In addition, comparison of our top candidate genes with two previously curated High and Low Activity gene lists highlight 24 overlapping genes, where Ndufa13, which encodes the supernumerary subunit A13 of complex I, was the only gene to be included in all three lists. Mitochondrial dysfunction has recently been implicated as both a cause and effect of anxiety-related disorders and thus should be further explored as a possible novel pharmaceutical treatment for anxiety disorders.

Evaluation and characterization of expression quantitative trait analysis methods in the Hybrid Rat Diversity Panel.

The Hybrid Rat Diversity Panel (HRDP) is a stable and well-characterized set of more than 90 inbred rat strains that can be leveraged for systems genetics approaches to understanding the genetic and genomic variation associated with complex disease. The HRDP exhibits substantial between-strain diversity while retaining substantial within-strain isogenicity, allowing for the precise mapping of genetic variation associated with complex phenotypes and providing statistical power to identify associated variants. In order to robustly identify associated genetic variants, it is important to account for the population structure induced by inbreeding. To this end, we investigate the performance of four plausible approaches towards modeling quantitative traits in the HRDP and quantify their operating characteristics. In particular, we investigate three approaches based on genome-wide mixed model analysis, and one approach based on ordinary least squares linear regression. Towards facilitating study planning and design, we conduct extensive simulations to investigate the power of genetic association analyses in the HRDP, and characterize the impressive attained power. In simulation of eQTL data in the HRDP, we find that a mixed model approach that leverages leave-one-chromosome-out kinship estimation attains the highest power while controlling type I error.

Loss of Frizzled 9 in Lung Cells Alters Epithelial Phenotype and Promotes Premalignant Lesion Development.

The transmembrane receptor Frizzled 9 (FZD9) is important for fetal neurologic and bone development through both canonical and non-canonical WNT/FZD signaling. In the adult lung, however, Fzd9 helps to maintain a normal epithelium by signaling through peroxisome proliferator activated receptor γ (PPARγ). The effect of FZD9 loss on normal lung epithelial cells and regulators of its expression in the lung are unknown. We knocked down FZD9 in human bronchial epithelial cell (HBEC) lines and found that downstream EMT targets and PPARγ activity are altered. We used a FZD9-/- mouse in the urethane lung adenocarcinoma model and found FZD9-/- adenomas had more proliferation, increased EMT signaling, decreased activation of PPARγ, increased expression of lung cancer associated genes, increased transformed growth, and increased potential for invasive behavior. We identified PPARγ as a transcriptional regulator of FZD9. We also demonstrated that extended cigarette smoke exposure in HBEC leads to decreased FZD9 expression, decreased activation of PPARγ, and increased transformed growth, and found that higher exposure to cigarette smoke in human lungs leads to decreased FZD9 expression. These results provide evidence for the role of FZD9 in lung epithelial maintenance and in smoking related malignant transformation. We identified the first transcriptional regulator of FZD9 in the lung and found FZD9 negative lesions are more dangerous. Loss of FZD9 creates a permissive environment for development of premalignant lung lesions, making it a potential target for intervention.

Changes in the Coexpression of Innate Immunity Genes During Persistent Islet Autoimmunity Are Associated With Progression of Islet Autoimmunity: Diabetes Autoimmunity Study in the Young (DAISY).

Longitudinal changes in gene expression during islet autoimmunity (IA) may provide insight into biological processes that explain progression to type 1 diabetes (T1D). We identified individuals from Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, autoantibodies present on two or more visits. Illumina's NovaSeq 6000 was used to quantify gene expression in whole blood. With linear mixed models we tested for changes in expression after IA that differed across individuals who progressed to T1D (progressors) (n = 25), reverted to an autoantibody-negative stage (reverters) (n = 47), or maintained IA positivity but did not develop T1D (maintainers) (n = 66). Weighted gene coexpression network analysis was used to identify coexpression modules. Gene Ontology pathway analysis of the top 150 differentially expressed genes (nominal P < 0.01) identified significantly enriched pathways including leukocyte activation involved in immune response, innate immune response, and regulation of immune response. We identified a module of 14 coexpressed genes with roles in the innate immunity. The hub gene, LTF, is known to have immunomodulatory properties. Another gene within the module, CAMP, is potentially relevant based on its role in promoting ß-cell survival in a murine model. Overall, results provide evidence of alterations in expression of innate immune genes prior to onset of T1D.

Relationship Between a Vitamin D Genetic Risk Score and Autoantibodies Among First-Degree Relatives of Probands With Rheumatoid Arthritis and Systemic Lupus Erythematosus.

Objective: Higher 25-hydroxyvitamin D (25(OH)D) levels have been associated with reduced risk for autoimmune diseases and are influenced by vitamin D metabolism genes. We estimated genetically-determined vitamin D levels by calculating a genetic risk score (GRS) and investigated whether the vitamin D GRS was associated with the presence of autoantibodies related to rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) in those at increased risk for developing RA and SLE, respectively. Methods: In this cross-sectional study, we selected autoantibody positive (aAb+) and autoantibody negative (aAb-) individuals from the Studies of the Etiologies of Rheumatoid Arthritis (SERA), a cohort study of first-degree relatives (FDRs) of individuals with RA (189 RA aAb+, 181 RA aAb-), and the Lupus Family Registry and Repository (LFRR), a cohort study of FDRs of individuals with SLE (157 SLE aAb+, 185 SLE aAb-). Five SNPs known to be associated with serum 25(OH)D levels were analyzed individually as well as in a GRS: rs4588 (GC), rs12785878 (NADSYN1), rs10741657 (CYP2R1), rs6538691 (AMDHD1), and rs8018720 (SEC23A). Results: Both cohorts had similar demographic characteristics, with significantly older and a higher proportion of males in the aAb+ FDRs. The vitamin D GRS was inversely associated with RA aAb+ (OR = 0.85, 95% CI = 0.74-0.99), suggesting a possible protective factor for RA aAb positivity in FDRs of RA probands. The vitamin D GRS was not associated with SLE aAb+ in the LFRR (OR = 1.09, 95% CI = 0.94-1.27). The SEC23A SNP was associated with RA aAb+ in SERA (OR = 0.65, 95% CI = 0.43-0.99); this SNP was not associated with SLE aAb+ in LFRR (OR = 1.41, 95% CI = 0.90 - 2.19). Conclusion: Genes associated with vitamin D levels may play a protective role in the development of RA aAbs in FDRs of RA probands, perhaps through affecting lifelong vitamin D status. The GRS and the SEC23A SNP may be of interest for future investigation in pre-clinical RA. In contrast, these results do not support a similar association in SLE FDRs, suggesting other mechanisms involved in the relationship between vitamin D and SLE aAbs not assessed in this study.

Epigenome-Wide Association Study of Infant Feeding and DNA Methylation in Infancy and Childhood in a Population at Increased Risk for Type 1 Diabetes.

We assessed associations between infant diet (e.g., breastfeeding and introduction to solid foods) and DNA methylation in infancy and childhood. We measured DNA methylation in peripheral blood collected in infancy (9-15 months of age) in 243 children; and in a subset of 50 children, we also measured methylation in childhood (6-9 years of age) to examine persistence, and at birth (in cord blood) to examine temporality. We performed multivariable linear regression of infant diet on the outcome of methylation using epigenome-wide and candidate site approaches. We identified six novel CpG sites associated with breastfeeding duration using an EWAS approach. One differentially methylated site presented directionally consistent associations with breastfeeding (cg00574958, CPT1A) in infancy and childhood but not at birth. Two differentially methylated sites in infancy (cg19693031, TXNIP; cg23307264, KHSRP) were associated with breastfeeding and were not present at birth; however, these associations did not persist into childhood. Associations between infant diet and methylation in infancy at three sites (cg22369607, AP001525.1; cg2409200, TBCD; cg27173510, PGBD5) were also present at birth, suggesting the influence of exposures other than infant diet. Infant diet exposures are associated with persistent methylation differences in CPT1A, which may be one mechanism behind infant diet's long-term health effects.

An effective processing pipeline for harmonizing DNA methylation data from Illumina's 450K and EPIC platforms for epidemiological studies.

OBJECTIVE: Illumina BeadChip arrays are commonly used to generate DNA methylation data for large epidemiological studies. Updates in technology over time create challenges for data harmonization within and between studies, many of which obtained data from the older 450K and newer EPIC platforms. The pre-processing pipeline for DNA methylation is not trivial, and influences the downstream analyses. Incorporating different platforms adds a new level of technical variability that has not yet been taken into account by recommended pipelines. Our study evaluated the performance of various tools on different versions of platform data harmonization at each step of pre-processing pipeline, including quality control (QC), normalization, batch effect adjustment, and genomic inflation. We illustrate our novel approach using 450K and EPIC data from the Diabetes Autoimmunity Study in the Young (DAISY) prospective cohort. RESULTS: We found normalization and probe filtering had the biggest effect on data harmonization. Employing a meta-analysis was an effective and easily executable method for accounting for platform variability. Correcting for genomic inflation also helped with harmonization. We present guidelines for studies seeking to harmonize data from the 450K and EPIC platforms, which includes the use of technical replicates for evaluating numerous pre-processing steps, and employing a meta-analysis.

A Mediation Approach to Discovering Causal Relationships between the Metabolome and DNA Methylation in Type 1 Diabetes.

Environmental factors including viruses, diet, and the metabolome have been linked with the appearance of islet autoimmunity (IA) that precedes development of type 1 diabetes (T1D). We measured global DNA methylation (DNAm) and untargeted metabolomics prior to IA and at the time of seroconversion to IA in 92 IA cases and 91 controls from the Diabetes Autoimmunity Study in the Young (DAISY). Causal mediation models were used to identify seven DNAm probe-metabolite pairs in which the metabolite measured at IA mediated the protective effect of the DNAm probe measured prior to IA against IA risk. These pairs included five DNAm probes mediated by histidine (a metabolite known to affect T1D risk), one probe (cg01604946) mediated by phostidyl choline p-32:0 or o-32:1, and one probe (cg00390143) mediated by sphingomyelin d34:2. The top 100 DNAm probes were over-represented in six reactome pathways at the FDR <0.1 level (q = 0.071), including transport of small molecules and inositol phosphate metabolism. While the causal pathways in our mediation models require further investigation to better understand the biological mechanisms, we identified seven methylation sites that may improve our understanding of epigenetic protection against T1D as mediated by the metabolome.

Severity of Idiopathic Scoliosis Is Associated with Differential Methylation: An Epigenome-Wide Association Study of Monozygotic Twins with Idiopathic Scoliosis.

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Inverse probability weighting is an effective method to address selection bias during the analysis of high dimensional data.

Omics studies frequently use samples collected during cohort studies. Conditioning on sample availability can cause selection bias if sample availability is nonrandom. Inverse probability weighting (IPW) is purported to reduce this bias. We evaluated IPW in an epigenome-wide analysis testing the association between DNA methylation (261,435 probes) and age in healthy adolescent subjects (n = 114). We simulated age and sex to be correlated with sample selection and then evaluated four conditions: complete population/no selection bias (all subjects), naïve selection bias (no adjustment), and IPW selection bias (selection bias with IPW adjustment). Assuming the complete population condition represented the "truth," we compared each condition to the complete population condition. Bias or difference in associations between age and methylation was reduced in the IPW condition versus the naïve condition. However, genomic inflation and type 1 error were higher in the IPW condition relative to the naïve condition. Postadjustment using bacon, type 1 error and inflation were similar across all conditions. Power was higher under the IPW condition compared with the naïve condition before and after inflation adjustment. IPW methods can reduce bias in genome-wide analyses. Genomic inflation is a potential concern that can be minimized using methods that adjust for inflation.

The oxylipin profile is associated with development of type 1 diabetes: the Diabetes Autoimmunity Study in the Young (DAISY).

AIMS/HYPOTHESIS: Oxylipins are lipid mediators derived from polyunsaturated fatty acids. Some oxylipins are proinflammatory (e.g. those derived from arachidonic acid [ARA]), others are pro-resolving of inflammation (e.g. those derived from α-linolenic acid [ALA], docosahexaenoic acid [DHA] and eicosapentaenoic acid [EPA]) and others may be both (e.g. those derived from linoleic acid [LA]). The goal of this study was to examine whether oxylipins are associated with incident type 1 diabetes. METHODS: We conducted a nested case-control analysis in the Diabetes Autoimmunity Study in the Young (DAISY), a prospective cohort study of children at risk of type 1 diabetes. Plasma levels of 14 ARA-derived oxylipins, ten LA-derived oxylipins, six ALA-derived oxylipins, four DHA-derived oxylipins and two EPA-related oxylipins were measured by ultra-HPLC-MS/MS at multiple timepoints related to autoantibody seroconversion in 72 type 1 diabetes cases and 71 control participants, which were frequency matched on age at autoantibody seroconversion (of the case), ethnicity and sample availability. Linear mixed models were used to obtain an age-adjusted mean of each oxylipin prior to type 1 diabetes. Age-adjusted mean oxylipins were tested for association with type 1 diabetes using logistic regression, adjusting for the high risk HLA genotype HLA-DR3/4,DQB1*0302. We also performed principal component analysis of the oxylipins and tested principal components (PCs) for association with type 1 diabetes. Finally, to investigate potential critical timepoints, we examined the association of oxylipins measured before and after autoantibody seroconversion (of the cases) using PCs of the oxylipins at those visits. RESULTS: The ARA-related oxylipin 5-HETE was associated with increased type 1 diabetes risk. Five LA-related oxylipins, two ALA-related oxylipins and one DHA-related oxylipin were associated with decreased type 1 diabetes risk. A profile of elevated LA- and ALA-related oxylipins (PC1) was associated with decreased type 1 diabetes risk (OR 0.61; 95% CI 0.40, 0.94). A profile of elevated ARA-related oxylipins (PC2) was associated with increased diabetes risk (OR 1.53; 95% CI 1.03, 2.29). A critical timepoint analysis showed type 1 diabetes was associated with a high ARA-related oxylipin profile at post-autoantibody-seroconversion but not pre-seroconversion. CONCLUSIONS/INTERPRETATION: The protective association of higher LA- and ALA-related oxylipins demonstrates the importance of both inflammation promotion and resolution in type 1 diabetes. Proinflammatory ARA-related oxylipins may play an important role once the autoimmune process has begun.

Phospholipid Levels at Seroconversion Are Associated With Resolution of Persistent Islet Autoimmunity: The Diabetes Autoimmunity Study in the Young.

Reversion of islet autoimmunity (IA) may point to mechanisms that prevent IA progression. We followed 199 individuals who developed IA during the Diabetes Autoimmunity Study in the Young. Untargeted metabolomics was performed in serum samples following IA. Cox proportional hazards models were used to test whether the metabolites (2,487) predicted IA reversion: two or more consecutive visits negative for all autoantibodies. We conducted a principal components analysis (PCA) of the top metabolites; |hazard ratio (HR) >1.25| and nominal P < 0.01. Phosphatidylcholine (16:0_18:1(9Z)) was the strongest individual metabolite (HR per 1 SD 2.16, false discovery rate (FDR)-adjusted P = 0.0037). Enrichment analysis identified four clusters (FDR P < 0.10) characterized by an overabundance of sphingomyelin (d40:0), phosphatidylcholine (16:0_18:1(9Z)), phosphatidylcholine (30:0), and l-decanoylcarnitine. Overall, 63 metabolites met the criteria for inclusion in the PCA. PC1 (HR 1.4, P < 0.0001), PC2 (HR 0.85, P = 0.0185), and PC4 (HR 1.28, P = 0.0103) were associated with IA reversion. Given the potential influence of diet on the metabolome, we investigated whether nutrients were correlated with PCs. We identified 20 nutrients that were correlated with the PCs (P < 0.05). Total sugar intake was the top nutrient. Overall, we identified an association between phosphatidylcholine, sphingomyelin, and carnitine levels and reversion of IA.

Androgen Receptor Regulates CD44 Expression in Bladder Cancer.

The androgen receptor (AR) is important in the development of both experimental and human bladder cancer. However, the role of AR in bladder cancer growth and progression is less clear, with literature indicating that more advanced stage and grade disease are associated with reduced AR expression. To determine the mechanisms underlying these relationships, we profiled AR-expressing human bladder cancer cells by AR chromatin immunoprecipitation sequencing and complementary transcriptomic approaches in response to in vitro stimulation by the synthetic androgen R1881. In vivo functional genomics consisting of pooled shRNA or pooled open reading frame libraries was employed to evaluate 97 genes that recapitulate the direction of expression associated with androgen stimulation. Interestingly, we identified CD44, the receptor for hyaluronic acid, a potent biomarker and driver of progressive disease in multiple tumor types, as significantly associated with androgen stimulation. CRISPR-based mutagenesis of androgen response elements associated with CD44 identified a novel silencer element leading to the direct transcriptional repression of CD44 expression. In human patients with bladder cancer, tumor AR and CD44 mRNA and protein expression were inversely correlated, suggesting a clinically relevant AR-CD44 axis. Collectively, our work describes a novel mechanism partly explaining the inverse relationship between AR and bladder cancer tumor progression and suggests that AR and CD44 expression may be useful for prognostication and therapeutic selection in primary bladder cancer. SIGNIFICANCE: This study describes novel AREs that suppress CD44 and an expected inverse correlation of AR-CD44 expression observed in human bladder tumors.

HDAC Inhibition Reverses Preexisting Diastolic Dysfunction and Blocks Covert Extracellular Matrix Remodeling.

BACKGROUND: Diastolic dysfunction (DD) is associated with the development of heart failure and contributes to the pathogenesis of other cardiac maladies, including atrial fibrillation. Inhibition of histone deacetylases (HDACs) has been shown to prevent DD by enhancing myofibril relaxation. We addressed the therapeutic potential of HDAC inhibition in a model of established DD with preserved ejection fraction. METHODS: Four weeks after uninephrectomy and implantation with deoxycorticosterone acetate pellets, when DD was clearly evident, 1 cohort of mice was administered the clinical-stage HDAC inhibitor ITF2357/Givinostat. Echocardiography, blood pressure measurements, and end point invasive hemodynamic analyses were performed. Myofibril mechanics and intact cardiomyocyte relaxation were assessed ex vivo. Cardiac fibrosis was evaluated by picrosirius red staining and second harmonic generation microscopy of left ventricle (LV) sections, RNA sequencing of LV mRNA, mass spectrometry-based evaluation of decellularized LV biopsies, and atomic force microscopy determination of LV stiffness. Mechanistic studies were performed with primary rat and human cardiac fibroblasts. RESULTS: HDAC inhibition normalized DD without lowering blood pressure in this model of systemic hypertension. In contrast to previous models, myofibril relaxation was unimpaired in uninephrectomy/deoxycorticosterone acetate mice. Furthermore, cardiac fibrosis was not evident in any mouse cohort on the basis of picrosirius red staining or second harmonic generation microscopy. However, mass spectrometry revealed induction in the expression of >100 extracellular matrix proteins in LVs of uninephrectomy/deoxycorticosterone acetate mice, which correlated with profound tissue stiffening based on atomic force microscopy. ITF2357/Givinostat treatment blocked extracellular matrix expansion and LV stiffening. The HDAC inhibitor was subsequently shown to suppress cardiac fibroblast activation, at least in part, by blunting recruitment of the profibrotic chromatin reader protein BRD4 (bromodomain-containing protein 4) to key gene regulatory elements. CONCLUSIONS: These findings demonstrate the potential of HDAC inhibition as a therapeutic intervention to reverse existing DD and establish blockade of extracellular matrix remodeling as a second mechanism by which HDAC inhibitors improve ventricular filling. Our data reveal the existence of pathophysiologically relevant covert or hidden cardiac fibrosis that is below the limit of detection of histochemical stains such as picrosirius red, highlighting the need to evaluate fibrosis of the heart using diverse methodologies.